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METHOD AND SYSTEM FOR OPTIMIZED 
PRE-S ATURATION IN MR WITH CORRECTED 
TRANSMITTER FREQUENCY OF PRE-PULSES 

BACKGROUND OF THE INVENTION 

[0001] The present invention relates to the field of Magnetic 
Resonance Imaging (MRI). More specifically, the invention relates to the field of 
pre saturation based MRI. 

[0002] Magnetic Resonance Imaging is a useful imaging tool for 
non-invasive imaging of internal organs. One of the important reasons for the rise in 
popularity of MRI is the outstanding clarity of the MRI images. Also, MRI uses radio 
frequency waves rather than X-rays. 

[0003] An MRI system works on the principle of magnetic resonance 
to obtain images of the human body. Human body tissues are composed of molecules 
like water and fat, which contain hydrogen atoms. The nuclei of hydrogen atoms have 
an associated magnetic moment, which is a measure of the net magnetic properties of 
the hydrogen nuclei. In a typical MRI system, electromagnets provide a strong 
magnetic field called the B.sub.O (B 0 ) magnetic field. The B 0 magnetic field interacts 
with the magnetic moments of hydrogen nuclei and causes a fraction of the nuclei to 
align with the B 0 magnetic field. The B 0 magnetic field also causes the magnetic 
moments of the aligned nuclei to precess around B 0 magnetic field direction at a 
frequency called frequency of precession (co). The frequency of precession depends 
on strength of the B 0 magnetic field and the gyromagnetic ratio of the nucleus (y), 
which is a characteristic property of the nucleus. Mathematically, co, B 0 and y are 
related as: 

<o = y*B 0 equation (1) 

[0004] The frequency of precession of the nucleus is called the 
Larmor frequency. The MRI system then applies radio frequency (RF) pulses at the 
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Larmor frequency of hydrogen nuclei. The frequency of the applied RF pulses is 
called the transmit frequency. When the precessing hydrogen nuclei are exposed to 
the RF pulses at their Larmor frequency, magnetic resonance occurs. When magnetic 
resonance occurs, the precessing hydrogen nuclei emit energy in the form of specific 
radio frequency signals. The radio frequency signals generated by the resonating 
hydrogen nuclei in water molecules (present in the scan volume) are processed by the 
MRI system to generate an MRI image. In many applications, the radio frequency 
signals generated by the resonating hydrogen nuclei in molecules other than water are 
undesired because they may reduce the quality of the MRI image. 

[0005] An MRI system images a volume of the human body called 
the scan volume. One of the ways to image a scan volume is by dividing the scan 
volume into a number of slices called scan slices. The scan slices are then imaged one 
by one. The images of these scan slices can then be combined to form an image of the 
scan volume. To successfully image a scan volume, the Bo magnetic field must be 
homogeneous across the scan volume. Bo magnetic field inhomogeneity across the 
scan volume can have an adverse impact on the quality of MRI images. Referring to 
equation (1), the transmit frequency (which is equal to the Larmor frequency) depends 
on the strength of the B 0 magnetic field. If the B 0 magnetic field is inhomogeneous, 
hydrogen nuclei at different locations of the scan volume will experience different Bo 
magnetic field strengths. This implies that the transmit frequency to be applied for 
magnetic resonance to occur will vary across the scan volume. A typical MRI system 
uses a single transmit frequency for the entire scan volume. Due to this, there will be 
some nuclei in the scan volume that will not undergo magnetic resonance. These 
nuclei will not produce a proper radio frequency signal, which leads to a poor quality 
MRI image. Hence, to obtain a good quality MRI image, a homogeneous B 0 magnetic 
field across the scan volume is desirable. 

[0006] A homogeneous Bo magnetic field is especially desirable in 
MRI systems that use pre-saturation. Pre-saturation is a method used for suppressing 
undesired signals that reduce image quality. In pre-saturation, frequency selective 
saturation pulses called RF pre-pulses are applied before the RF pulses. The use of 
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frequency selective RF pre-pulses suppresses the undesired signals. An example of 
such RF pre-pulses is fat saturation. Fat saturation is a technique that selectively 
suppresses undesired signals from hydrogen nuclei present in fat molecules. Selective 
suppression is achieved by applying specific fat saturation RF pre-pulses prior to the 
RF pulses. 

[0007] Any inhomogeneity in the B 0 magnetic field can lead to 
wrong determination of frequency of RF pre-pulses with respect to the RF pulse used 
to excite hydrogen nuclei present in water molecules. The use of RF pre-pulses at a 
wrong frequency substantially affects the image quality. For example in fat saturation, 
the RF pre-pulses at a wrong frequency will not completely suppress the signals from 
hydrogen nuclei present in fat molecules. In some cases, the fat saturation pulses of a 
wrong frequency may suppress signals from hydrogen nuclei present in water 
molecules, thereby reducing the image quality. Fig. 1 shows the effect of B 0 magnetic 
field inhomogeneity on an MRI image. The MRI image shows poor fat saturation and 
unwanted suppression of signals from hydrogen nuclei present in water molecules. 

[0008] Hence, there exists a need for a method to reduce the effect of 
B 0 magnetic field inhomogeneity in pre-saturation based MRI experiments. 

BRIEF SUMMARY OF THE INVENTION 

[0009] In accordance with one aspect of the present invention, a 
method is provided for optimizing pre-saturation in magnetic resonance imaging. The 
method first generates a B 0 field map of each scan slice of a scan volume. Thereafter, 
a first frequency of RF pre-pulses is obtained by a standard procedure. The Bo field 
map of each scan slice is then used to calculate the median value of the Bo magnetic 
field over each scan slice as well as the percentage of positive and negative scan slice 
pixels in each scan slice. In case the percentage of either the positive scan slice pixels 
or the negative scan slice pixels in each scan slice is greater than a predefined 
threshold value, a second frequency of RF pre-pulses is calculated for each scan slice. 
The second frequency of RF pre-pulses for a scan slice is calculated by adding the 
median value of the B 0 magnetic field over the scan slice to the first frequency of RF 
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pre-pulses. However, if the percentage value of positive scan slice pixels or the 
percentage value of negative scan slice pixels of a scan slice does not exceed the 
predetermined threshold value, a user of the MRI system is advised to improve 
shimming. After shimming has been improved by the user, the above-mentioned 
steps are repeated for each scan slice, until the percentage of either the positive scan 
slice pixels or the negative scan slice pixels exceeds the predefined threshold value. 
Thereafter, MRI images of the scan slices are obtained using the RF pre-pulses at the 
second frequency for each scan slice. 

[0010] In accordance with another aspect of the present invention, an 
MRI system is provided for generating MRI images of a scan volume. The MRI 
system comprises a magnet for producing a high intensity magnetic field called the Bo 
magnetic field. The MRI system further comprises a set of shimming coils for 
improving Bo magnetic field homogeneity and a set of gradient coils for producing 
gradient magnetic fields superposed on the Bo magnetic field. The gradient magnetic 
fields are used by the MRI system to select a specific region of the scan volume for 
magnetic resonance imaging. A database is provided, which stores the B 0 magnetic 
field distribution over each scan slice in the form of B 0 field maps. A processing 
module is provided for calculating a first frequency of RF pre-pulses for the scan 
slice. The processing module further calculates the median value of the Bo magnetic 
field over each scan slice. The processing module further calculates a second 
frequency of RF pre-pulses for each scan slice by adding the median value of the Bo 
magnetic field over the scan slice to the first frequency of RF pre-pulses. A 
transmitter is provided, which generates RF pre-pulses at the calculated second 
frequency for pre-saturation and RF pulses for causing magnetic resonance of the 
hydrogen nuclei in the scan volume. A detector is provided for detecting radio 
frequency signals generated due to magnetic resonance. The processing module 
processes the radio frequency signals generated due to magnetic resonance to obtain 
MRI images. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[001 1] Various embodiments of the present invention will hereinafter 
be described in conjunction with the appended drawings provided to illustrate and not 
to limit the present invention, wherein like designations denote like elements, and in 
which: 

[0012] Fig. 1 is an MRI image showing poor fat saturation and 
unwanted suppression of signals from hydrogen nuclei present in water molecules due 
to B 0 magnetic field inhomogeneity; 

[0013] Figs. 2 A and 2B illustrate a flowchart depicting the steps 
involved in the method of optimizing pre-saturation in magnetic resonance imaging, 
in accordance with one embodiment of the present invention; 

[0014] Fig. 3 is an MRI image obtained from an MRI system using 
the method of the present invention; 

[0015] Fig. 4 is a block diagram showing the general working 
environment of an MRI system for pre-saturation based MRI in accordance with one 
embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

[0016] The present invention relates to a method, system and 
computer program product for optimizing pre-saturation in magnetic resonance 
imaging. The optimization is done by setting a correct transmitter frequency of RF 
pre-pulses for all scan slices. 

[0017] In one embodiment, the present invention is envisioned to be 
operating in conjunction with an MRI system using the technique of fat saturation. 
Fat saturation is a technique used to reduce interference from magnetic resonance 
signals generated by hydrogen nuclei present in fat molecules. 
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[0018] It will be apparent to a person skilled in the art that the 
present invention can be adapted to operate with MRI systems that use pre-saturation 
techniques other than fat saturation. For example, the present invention can find 
application in magnetization transfer for suppression of signals from macromolecules 
(such as proteins) and in spectroscopy experiments for suppression of signals from 
water molecules. 

[0019] Figs. 2 A and 2B illustrate a flowchart depicting the steps 
involved in the method of optimizing pre-saturation in magnetic resonance imaging, 
in accordance with one embodiment of the present invention. 

[0020] At step 202, distribution of the B 0 magnetic field over each 
scan slice of the scan volume is obtained. There exist several methods to obtain the 
B 0 magnetic field distribution. Some examples of such methods include dual gradient 
echo method, dual spin echo method and dual echo spiral scan. The Bo magnetic field 
distribution is obtained as a Bo field map of each scan slice. The Bo field map of a 
scan slice is a map of the Bo magnetic field over the scan slice. Value of a pixel of the 
Bo field map of a scan slice is the frequency offset of the B 0 magnetic field over the 
pixel from the average B 0 magnetic field. In one embodiment of the present 
invention, the B 0 field maps are stored in a database. 

[0021] At step 204, a first frequency for the RF pre-pulses is obtained 
by a standard procedure. The procedure involves calculating a transmit frequency for 
the RF pulses. This is done by an experiment that calculates the frequency at which 
the maximum hydrogen nuclei present in water molecules resonate. Once the transmit 
frequency has been calculated, the first frequency of the RF pre-pulse is calculated 
based on the difference between the frequencies at which fat and water molecules 
resonate at the average B 0 magnetic field. For example, at an average Bo magnetic 
field of 1 .5 Tesla, this frequency difference is 220 Hz. Thus, if the transmit frequency 
is set to 63,584,500 Hz then the RF pre-pulse frequency will be set to 63,584,280 Hz. 
The frequency 63,584,280 Hz will be the first frequency of RF pre-pulses. 
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[0022] The type and frequency of RF pre-pulses depend on the type 
of molecules whose magnetic resonance signals are to be suppressed. In one 
embodiment of the invention, the RF pre-pulses are used to suppress magnetic 
resonance signals from fat molecules. 

[0023] At step 206, a median value of the B 0 magnetic field over 
each scan slice is calculated. Calculation of the median value of the B 0 magnetic field 
over a scan slice is done by calculating the median of the Bo map pixel values of the 
scan slice. In one embodiment, the median value of the Bo magnetic field over each 
scan slice is calculated using the B 0 field maps stored in a database. 

[0024] At step 208, the percentage of positive scan slice pixels and 
negative scan slice pixels in each scan slice is calculated. The calculation is done for 
a scan slice in the following manner. First, the total number of positive and negative 
scan slice pixels in the scan slice is calculated. A positive scan slice pixel is defined 
as a pixel that has a positive or zero value on the B 0 field map of the scan slice. 
Similarly, a negative scan slice pixel is defined as a pixel that has a negative value on 
the Bo field map of the scan slice. Thereafter, the total number of pixels in the scan 
slice is calculated. 

[0025] The percentage of positive scan slice pixels in a scan slice is 
obtained by dividing the number of positive scan slice pixels in the scan slice by the 
total number of pixels in the scan slice. Similarly, the percentage of negative scan 
slice pixels in a scan slice is obtained by dividing the number of negative scan slice 
pixels in the scan slice by the total number of pixels in the scan slice. In one 
embodiment of the present invention, the percentage of positive scan slice pixels and 
negative scan slice pixels in each scan slice is calculated from the Bo magnetic field 
maps stored in a database. 

[0026] At step 210, a check is made for each scan slice to determine 
if either the percentage value of positive scan slice pixels or the percentage value of 
negative scan slice pixels of the scan slice, exceed a predetermined threshold value. 
In one embodiment of the present invention, the threshold value is 80%. 
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[0027] At step 210, if for at least one scan slice, the percentage value 
of positive scan slice pixels or the percentage value of negative scan slice pixels of the 
scan slice does not exceed the predetermined threshold value, step 212 is performed. 

[0028] At step 212, the user of the MRI system is advised to improve 
the shimming. It would be evident to one skilled in the art that there are many ways of 
improving the shimming. For example, if the MRI system is using a shim protocol 
that corrects only linear gradients in the Bo magnetic field, the user may be advised to 
improve shimming by using a higher order shim protocol. The user can then improve 
the shimming at step 213 by following the advice. 

[0029] After the shimming has been improved at step 213, the 
method returns to step 202. 

[0030] Referring back to step 210, if for each scan slice, the 
percentage value of positive scan slice pixels or the percentage value of negative scan 
slice pixels of the scan slice exceeds a predetermined threshold value, step 214 is 
performed. 

[0031] At step 214, a second frequency of RF pre-pulses is calculated 
for each scan slice. The calculation of the second frequency of RF pre-pulses for a 
scan slice is done by using the median value of the B 0 magnetic field over the scan 
slice and the first frequency of RF pre-pulses. In one embodiment, the second 
frequency of RF pre-pulses for a scan slice is calculated by adding the median value of 
the B 0 magnetic field over the scan slice to the first frequency of RF pre-pulses. 

[0032] At step 216, the scan slices are imaged. A scan slice is 
imaged by first applying RF pre-pulses at the second frequency calculated for the scan 
slice. Thereafter, RF pulses at transmit frequency are applied to the scan slice. 
Application of the RF pulses causes magnetic resonance of hydrogen nuclei present in 
the scan slice. When magnetic resonance occurs, the hydrogen nuclei generate radio 
frequency signals. The magnetic resonance signals are measured. The radio 
frequency signals are processed to obtain an MRI image of the scan slice. The use of 
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RF pre-pulses at the second frequency suppresses unwanted radio frequency signals 
from molecules other than water molecules. In one embodiment of the present 
invention, the RF pre-pulses are used for fat saturation. In such an embodiment, the 
RF pre-pulses suppress unwanted magnetic resonance signals from hydrogen nuclei 
present in fat molecules. 

[0033] At step 218, the MRI images obtained at step 216 are stored 
in a database. The MRI images stored at step 218 can then be displayed on a display 
device. 

[0034] Fig. 3 is an MRI image obtained from an MRI system using 
the method of the present invention. The MRI image of Fig. 3 shows improved fat 
saturation and reduced suppression of signals from hydrogen nuclei present in water 
molecules as compared to the MRI image in Fig. 1 . 

[0035] Fig. 4 is a block diagram showing the general working 
environment of an MRI system 400 for pre-saturation based MRI in accordance with 
one embodiment of the present invention. 

[0036] MRI system 400 further comprises a processing module 402. 
Processing module 402 can be in the form of a general-purpose computer, a 
programmed microprocessor, a micro-controller, a peripheral integrated circuit 
element, and other devices or arrangements of devices. MRI system 400 further 
comprises a control panel 404 connected to processing module 402. A user can 
control the MRI system by entering commands from control panel 404. 

[0037] MRI system 400 further comprises a database 406 and a 
display device 408 connected to processing module 402. Typical examples of display 
device 408 include liquid crystal displays (LCDs), cathode ray tubes (CRTs), light- 
emitting diode (LED) displays and television screens. 

[0038] MRI system 400 further comprises a strong, polarizing 
magnet 410 and a set of shim coils 412, both connected to processing module 402. 
Magnet 410 generates a strong magnetic field called B 0 magnetic field over a scan 
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volume. Shim coils 412 generate an additional magnetic field that reduces 
inhomogeneities in the magnetic field Bo over the scan volume. 

[0039] MRI system 400 further comprises a magnetic field detector 
414 connected to processing module 402. Magnetic field detector 414 measures the 
Bo magnetic field distribution over all scan slices of the scan volume and sends the Bo 
magnetic field distribution data to processing module 402. Processing module 402 
then generates a Bo field map of all scan slices from the Bo magnetic field distribution 
data. The Bo field map of a scan slice is a map of the Bo magnetic field over the scan 
slice. Thereafter, processing module 402 stores the Bo field maps in database 406. 
Further, processing module 402 calculates the median of the Bo magnetic field over 
each scan slice by calculating the median of pixel values of the Bo field map. 
Thereafter, processing module 402 stores the calculated median values for each scan 
slice in database 406. Processing module 402 further calculates a second frequency of 
RF pre-pulses for each scan slice by correcting a first frequency of RF pre-pulses for 
the scan slice. (The first frequency of RF pre-pulses is obtained by a standard 
procedure). In one embodiment, the correction is done by adding the median value of 
the B 0 magnetic field over the scan slice to the first frequency of RF pre-pulses. 

[0040] MRI system 400 further comprises a set of gradient coils 416 
connected to processing module 402. Gradient coils 416 are current carrying coils 
that produce gradient magnetic fields. These gradient magnetic fields are superposed 
on the Bo magnetic field by MRI system 400 to select a specific region of the scan 
volume for magnetic resonance imaging. In one embodiment, gradient coils 416 are 
three paired orthogonal current-carrying coils. 

[0041] MRI system 400 further comprises a transmitter 418 
connected to processing module 402. Transmitter 418 comprises a coil that generates 
RF pulses and RF pre-pulses. 

[0042] The RF pulses (at the transmit frequency) generated by 
transmitter 418 cause magnetic resonance in hydrogen nuclei. Magnetic resonance 
causes hydrogen nuclei to generate radio frequency signals. The radio frequency 
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signals generated by the hydrogen nuclei are detected by a receiver 420, connected to 
processing module 402. Receiver 420 comprises a coil or antenna, positioned within 
the scan volume and connected to the receiver circuitry to detect the radio frequency 
signals generated due to magnetic resonance. Receiver coils are of various types, for 
example, solenoidal, planar, volume, quadrature and phased array coils. Receiver 420 
sends the detected radio frequency signals to processing module 402. Processing 
module 402 then processes the detected radio frequency signals and generates MRI 
images. 

[0043] The MRI images generated by processing module 402 are 
stored in database 406. The images stored in database 406 are displayed on displaying 
device 408. 

[0044] The system, as described in the present invention or any of its 
components may be embodied in the form of a processing machine. Typical examples 
of a processing machine include a general-purpose computer, a programmed 
microprocessor, a micro-controller, a peripheral integrated circuit element, and other 
devices or arrangements of devices, which are capable of implementing the steps that 
constitute the method of the present invention. 

[0045] The processing machine executes a set of instructions in the 
form of a computer program code. The computer program code can be stored in one 
or more storage units in order to process input data. The storage units may be in the 
form of a database or a physical memory element present in the processing machine. 

[0046] The computer program code may include various instructions 
that instruct the processing machine to perform specific tasks such as the steps that 
constitute the method of the present invention. The computer program code may be in 
various forms, such as system software or application software. Further, the computer 
program code can be in the form of a collection of separate programs, a program 
module with a larger program or a portion of a program module. The computer 
program code can be written in any programming language such as C, C++, C#, Java 
etc. The processing of input data by the processing machine may be in response to 
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user commands, or in response to results of previous processing or in response to a 
request made by another processing machine. 

[0047] It will be apparent to a person skilled in the art that it is not 
necessary for the various processing machines and/or storage units to be physically 
located in the same geographical location. The processing machines and/or storage 
units can be located in geographically distinct locations and connected to each other to 
enable communication. Various communication technologies can be used to enable 
communication between the processing machines and/or storage units. Such 
communication technologies can use various protocols such as TCP/IP, UDP, ATM or 
OSI. Such technologies include connection of the processing machines and/or storage 
units, in the form of a network. The network can be an intranet, an extranet, the 
Internet or any client server models that enable communication. 

[0048] In the system and method of the present invention, a variety of 
"user interfaces" can be utilized to allow a user to interface with the processing 
machine or machines that are used to implement the present invention. The user 
interface is used by the processing machine to interact with a user in order to convey 
or receive information. The user interface can be any hardware, software, or a 
combination of hardware and software used by the processing machine that allows a 
user to interact with the processing machine. The user interface can be in the form of 
a dialogue screen and can include various associated devices to enable communication 
between a user and a processing machine. It is contemplated that the user interface 
can interact with another processing machine rather than a human user. Further, it is 
also contemplated that the user interface can interact partially with other processing 
machines while also interacting partially with the human user. 

[0049] While the various embodiments of the present invention have 
been illustrated and described, it will be clear that the present invention is not limited 
to these embodiments only. Numerous modifications, changes, variations, 
substitutions and equivalents will be apparent to those skilled in the art without 
departing from the spirit and scope of the present invention as described in the claims. 
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